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ABSTRACT 

Galaxies are missing most of their baryons, and many models predict these 
baryons lie in a hot halo around galaxies. We establish observationally motivated 
constraints on the mass and radii of these haloes using a variety of independent 
arguments. First, the observed dispersion measure of pulsars in the Large Mag- 
ellanic Cloud allows us to constrain the hot halo around the Milky Way: if it 
obeys the standard NFW profile, it must contain less than 4-5% of the missing 
baryons from the Galaxy. This is similar to other upper limits on the Galactic 
hot halo, such as the soft X-ray background and the pressure around high veloc- 
ity clouds. Second, we note that the X-ray surface brightness of hot haloes with 
NFW profiles around large isolated galaxies is high enough that such emission 
should be observed, unless their haloes contain less than 10-25% of their missing 
baryons. Third, we place constraints on the column density of hot haloes using 
nondetections of OVII absorption along AGN sightlines: in general they must 
contain less than 70% of the missing baryons or extend to no more than 40 kpc. 
Flattening the density profile of galactic hot haloes weakens the surface bright- 
ness constraint so that a typical L* galaxy may hold half its missing baryons 
in its halo, but the OVII constraint remains unchanged, and around the Milky 
Way a flattened profile may only hold 6 — 13% of the missing baryons from the 
Galaxy (2 — 4 x 1O 1O M ). We also show that AGN and supernovae at low to 
moderate redshift - the theoretical sources of winds responsible for driving out 
the missing baryons - do not produce the expected correlations with the baryonic 
Tully-Fisher relationship and so are insufficient to explain the missing baryons 
from galaxies. We conclude that most of missing baryons from galaxies do not 
lie in hot haloes around the galaxies, and that the missing baryons never fell 
into the potential wells of protogalaxies in the first place. They may have been 
expelled from the galaxies as part of the process of galaxy formation. 
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1. Introduction 

The so-called "missing baryon problem" actually refers to two separate, but related, 
issues. First is the realization that the mean density of all detected baryons in the local 
universe accounts for less than half of the cosmological baryon density. The universal baryon 
fraction is known precisely from, among other sources, the Wilkinson Microwave Anisotropy 
Probe 5-year data (Dunkley et al. 2009): f b = Vt b /Vt m = 0.171 ±0.006 and Q b h 2 = 0.0227 ± 
0.0006. But studies that attempt to count baryons in the low-redshift universe (Fukugita, 
Hogan, and Peebles 1998, Fukugita and Peebles 2004, Nicastro et al. 2005a) find that only 
about a tenth of this figure is observed in stars and cool gas in galaxies, while another third 
of Q b is divided between the hot plasma in the intracluster medium of galaxy clusters, and 
the cool intercluster gas detected in Lya absorption lines (for more on the latter, see e.g. 
Penton, Stocke, and Shull 2004). 

This result introduces the second missing baryon problem, which is the realization that 
most galaxies are severely baryon-depleted relative to the cosmological fraction (e.g. Bell 
et al. 2003). Based on cosmological simulations, the rest of the baryon budget is thought 
to reside in the "warm-hot intergalactic medium" (WHIM), at temperatures between 10 5 
and 10 7 K (Cen and Ostriker 1999, Dave et al. 2001). This is thought to include most 
of the baryons associated with galaxies, since the optically-luminous portion of galaxies is 
severely baryon-depleted relative to the universal fraction (Dai et al. 2010). Observational 
evidence for the existence of the WHIM has subsequently begun to accrue, primarily through 
detections with ultraviolet and X-ray telescopes of highly-ionized Oxygen absorption lines 
(for a recent review, see Bregman 2007). 

These observations do not yet constrain the nature and distribution of the WHIM. 
Essentially all models agree that the WHIM exists as a cosmic filamentary web with some 
material in hot haloes around galaxies, but the mass fraction as a function of WHIM density 
is still unknown. Some numerical models (e.g Dave et al. 2009) find a diffuse WHIM at 
z = with most of its baryonic content lying outside of galactic haloes. However, many 
simulations either assume (e.g. Bower et al. 2006, Croton et al. 2006) or self-consistently 
derive (e.g. Cen and Ostriker 2006, Tang et al. 2009, Kim, Wise, and Abel 2009) the result 
that cosmologically significant reservoirs of baryons are embedded in hot haloes around 
massive galaxies (kT ~ a few hundred eV). These hot haloes have also been modeled in 
separate theoretical calculations (Fukugita and Peebles 2006, Dekel and Birnboim 2006, 
Kaufmann et al. 2009) and have important consequences in the galaxy formation process 
such as differentiating the "blue cloud" from the "red sequence" (Dekel and Birnboim 2006, 
Bouche et al. 2009). Observationally, the evidence for hot haloes is still unclear: gas at 
these temperatures has been observed in X-rays around disk galaxies out to a few kpc 
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beyond the disk (Strickland et al. 2004, Li et al. 2006, Li, Wang, and Hameed 2007), but 
this gas seems to result from starburst activity in the host galaxy instead of emerging as 
a byproduct of galaxy formation (Rasmussen 2009), and its inferred mass falls significantly 
short of theoretical predictions (Tullmann et al. 2006). 

In this paper, we examine the implications of this paradigm in more detail. If most 
of the baryons associated with galaxies reside in hot haloes around the optically luminous 
part, there are several predictions that can be tested with existing observations. We seek to 
answer two questions. First, how much baryonic mass resides in these haloes? And, what 
can these haloes and the second missing baryon problem tell us about the process of galaxy 
formation? 

To answer these questions, we present several independent lines of argument. First, we 
consider observational constraints on the density of hot gas around the Milky Way using the 
dispersion measure of pulsars in the Large Magellanic Cloud, as well as the ambient pressure 
around high velocity clouds, and the galactic soft X-ray background. We then extend this 
argument to other galaxies. We constrain the characteristic density and radius of hot haloes 
using existing X-ray observations of absorption along quasar sightlines and emission from 
galactic haloes. We also discuss the energetics of driving a galactic-scale wind, arguing that 
existing mechanisms for expelling most of the missing baryons seem incomplete. Finally, 
we discuss flattened density profiles and other prospects for resolving the second missing 
baryon problem. Before engaging these arguments, however, we will begin by describing the 
assumptions we use to model hot haloes, which are very similar to those of Fukugita and 
Peebles (2006). 



2. On Milky Way Parameters 

Much of the following discussion directly relates to the Milky Way or uses scaling re- 
lations to extrapolate from the Galaxy to other galaxies. Therefore we begin by discussing 
the parameters and scaling relations used throughout this paper. 

Recent estimates for the virial mass of the Milky Way vary significantly. The consensus 
figure is roughly 2.5 x 1O 12 M (Sakamoto, Chiba, and Beers 2003, Bellazzini 2004, Loeb et 
al. 2005, Li and White 2008) , although a recent estimate was as low as 1.0±° o i x 1O 12 M 
(Xue et al. 2008) while another (Reid et al. 2009) revised upwards the circular velocity by 
about 15 km s _1 , which should increase the mass estimate by more than 10%. Since larger 
masses for the Galaxy increase the amount of missing baryons, we will adopt a conservative 
virial mass estimate of 2.0 x lO^M^. 
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Regarding the virial radius of the Milky Way, there is scatter around a central value. 
Shattow and Loeb (2009) use i? v ; r = 277 for a solar distance from the Milky Way center of 
d & = 8.0 kpc, while the earlier work of Loeb et al. (2005) uses a virial radius of 207 kpc. 
Klypin, Zhao, and Somerville (2002) favor a virial radius of 258 kpc. We adopt an estimate 
for the virial radius of 250 kpc. 

Another relevant Galactic parameter is amount of missing baryons. There is a sur- 
prisingly tight correlation between the total mass of a galaxy and its baryon fraction (Mc- 
Gaugh 2000, Dai et al. 2010), which we parametrize as f b = 0.04(M/2 x lO^M©) 1 / 2 for 
M < 6 x 1O 12 M . The relation flattens significantly for galaxy groups and clusters, but those 
are outside the scope of this work. This implies a baryonic mass of 8 x 1O 1O M for the 
Galaxy, which broadly comports with the estimate of a stellar mass of 5 x 10 10 M Q for the 
Galaxy (Binney and Tremaine 2008), with an additional 1 — 2 x 10 10 M Q in gas and dust. We 
assume the total mass is comprised of three components: the dark matter halo, the observ- 
able baryons in the galaxy, and the missing baryons in a hot halo. In this model, the mass 
of the missing baryons can be related to the observed mass M of the other two components: 



For the Milky Way, the missing baryonic mass is therefore 3 x 1O 11 M . As a simple first- 
order model, we assume some fraction of the missing mass is distributed in a hot virialized 
halo around the disk. We can model the gaseous halo by assuming the Galactic dark matter 
halo obeys a Navarro, Frenk, and White (NFW) profile (Navarro, Frenk, and White 1997)). 
The missing baryons would follow the gravitational potential set by the dark matter, and 
therefore also obey the same density profile, except that the missing baryons do not extend 
into the disk, so we occasionally truncate the gaseous halo at an inner radius r 1 (this does 
not have a significant effect on our analysis). We examine many values for r\ and for the 
NFW concentration parameter C = ro/R vir , but our nominal choice will be r\ = 10 kpc and 
C = 12 (so r = 20.8 kpc). 

Finally, the gaseous halo is assumed to be approximately isothermal at the virial tem- 
perature kT vir ~ 300 eV. This is also the temperature of the hot gas observed around NGC 
891 (Tiillmann et al. 2006), a near 'twin' of the Milky Way (van der Kruit 1984), and is 
similar to the 250 eV gas surrounding the Galactic disk observed by Snowden et al. (1998). 
We expect the gas to have Z ~ 0.2Z Q (using Anders and Grevesse 1989 abundances) follow- 
ing Cen and Ostriker (2006) and Rasumssen et al. (2009), and since similar metallicity is 
found in the hot intracluster medium even towards the virial radius of distant galaxy clusters 
(Maughan et al. 2008, Anderson et al. 2009). In computing the electron density of the halo, 
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we take ji = 0.6 and \i e = 1.2. 



3. Constraints on the Milky Way Halo 



3.1. LMC Pulsar Dispersion Measure 



The strongest constraint we can derive for the hot halo around the Galaxy comes from 
the dispersion measure (DM) of pulsars in the Large Magellanic Cloud (LMC). The DM is 
defined as 



where d is the distance to the pulsar and n(l) is the free electron density along the line of 
sight. The Galactic free electrons (primarily from the so-called "Warm Ionized Medium", or 
"WIM" ) have been modeled as an exponential disk by, e.g, Gaensler et al. (2008). Integrating 
over this disk out to the LMC (and assuming a distance of 50 kpc from the Sun), we estimate 
the Galactic WIM contributes about 47 cm" 3 pc to the DM of pulsars in the LMC. In a 
survey of the Magellanic Clouds, Manchester et al. (2006) measured the DM of 11 pulsars in 
the direction of the LMC. Three pulsars have DMs below 47 cm" 3 pc (the values are 26,28, 
and 45 cm" 3 pc), suggesting these sources are located within the Galaxy and randomly 
superimposed in front of the LMC. The next-lowest DMs for LMC pulsars are 66, 69, and 
73 cm" 3 pc, which leads to a conservative estimate of ~ 70 cm" 3 pc for the minimum DM 
introduced by the free electrons towards the LMC. This is also consistent with the lowest 
DMs for pulsars in the same survey towards the Small Magellanic Cloud (SMC): 70 and 76 



Subtracting 47 cm" 3 pc of Galactic electrons leaves 23 cm" 3 pc of electrons from other 
sources. The LMC has an ionized halo that may be responsible for much of this excess DM. 
One can assess the amount of ionized gas from the high-ionization absorption lines, such as 
OVI and CIV. The column densities of these ions in the LMC are comparable to the Milky 
Way (Howk et al. 2002). Considering that the metallicity is lower in the LMC, the inferred 
electron column is likely to be greater than that of the Milky Way. This implies that only a 
fraction of the excess DM is due to an ambient medium along the line of sight toward the 
LMC. However, since we do not know precisely the position of the LMC pulsars within its 
halo, in order to keep our estimate as conservative as possible we neglect the LMC electrons 
and assign all 23 cm -3 pc of excess DM to the hot halo of the Milky Way. 
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With an LMC distance of 50 kpc, 23 cm -3 pc corresponds to an average electron density 
between the Sun and the LMC of < n >m 5 x 10" 4 cm" 3 . If the entire halo had this mean 
density out to the virial radius, it would contain 8 x 10 11 M Q - all of the missing baryons. 
This is an unphysical profile, however. We instead consider the standard assumption of 
an NFW profile for the missing baryons, with a concentration of 12 (Klypin, Zhao, and 
Somerville 2002), and find a halo obeying this profile and the DM constraint would have a 
mass of 1.2 x 1O 1O M , which is only 4% of the missing baryons from the Galaxy. However, 
since the Galactic electrons have a scale length of 1.8 kpc (Gaensler et al. 2008), we might 
reasonably argue that baryons within a few scale lengths are Galactic, not missing baryons. 
We therefore also consider an NFW profile with an inner truncation radius of 10 kpc, which 
is five scale lengths and therefore a very conservative estimate. Such a profile would have a 
mass of 1.5 x 10 10 M Q , or just 5% of the missing baryons from the Galaxy. 

As a consistency check, we compute the cooling time of this halo, using our halo temper- 
ature and metallicity with the cooling function of Sutherland and Dopita (1993). Our model 
assumes a stable halo, but we find the halo is probably not quite stable against cooling over 
a Hubble time. For n = 5 x 10" 4 cm" 3 (the upper limit on the average gas density out to 
the LMC), the cooling time is r ^ 10 Gyr. For the above NFW profile, then, 6.4 x 1O 9 M 
(40% of the total hot halo gas mass) will cool within the age of the Galaxy and re-accrete 
onto the disk. To maintain self-consistency, the constraints above must be treated as upper 
limits on an even more rarefied halo, or there must be continuous re-heating of the halo gas 
by a wind or some other mechanism. 



3.2. Other Constraints 

We can compare the constraints on the Galactic hot halo from LMC pulsars to other 
constraints derived using a variety of methods. For reference, the DM upper limit on the 
electron density at r = 50 kpc is n e < 7.7 x 10" 5 cm -3 for an NFW halo, or n e < 1.0 x 10" 4 
cm -3 if we truncate the profile at r in = 10 kpc. 

Our first comparison is to the result of Stanimirovic et al. (2002), who used 21 cm 
observations of the Magellanic Streams to measure the radius and velocity width of clumps 
of neutral gas in the streams. From these measurements, they derive the internal pressure 
of the clumps, and equate this to the gas pressure of the hot halo (making the assumption 
that the clumps are confined entirely by gas pressure, and neglecting magnetic and turbulent 
pressure; this is therefore probably an upper limit on the gas pressure). They assume a halo 
temperature of 1 x 10 6 K; for our higher temperature of 3.5 x 10 6 K, the inferred upper limit 
on the hot halo electron density is n(45 kpc)< 9 x 10~ 5 cm" 3 - consistent with the above 
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constraint. 

We also consider high velocity clouds, using UV measurements of the temperature and 
density of numerous partially- ionized metals and inferred gas pressures (Fox et al. 2005). 
This analysis suffers from a size-distance degeneracy, however, and so Fox et al. perform 
their analysis for several potential cloud distances. The lower distance (10 kpc) yields a 
lower upper limit, and the higher distance (100 kpc) yields a higher upper limit, but they 
seem to prefer a distance of 50 kpc based on a tentative association of some of these clouds 
with the Magellanic Stream. Using this distance and converting to our halo temperature, 
the upper limits on the electron density at this radius are approximately n e = 3 — 6 x 10~ 5 
cm' 3 . 

We also estimate density using the emission measure of the hot halo, as measured 
by Kuntz and Snowden (2000). They decompose the ROSAT soft X-ray background into 
components, tentatively identifying the harder emission from an unidentified component as 
the hot halo of missing baryons. This emission has a temperature of about 3 x 10 6 K, 
close to our estimate of the hot halo temperature, and it appears to have a fairly uniform 
distribution across the sky. There is also a softer background component with a temperature 
of T ~ 1 x 10 6 K, but this temperature is too cool to remain in a hot halo out to the virial 
radius, and its spatial distribution is not isotropic: it follows the Galactic disk. 

The hotter component of the Galactic X-ray background has a best-fit emission measure 
of 2.1 x 10~ 3 cm -6 pc, for an assumed solar-metallicity gas. This corresponds to an NFW 
hot halo of total mass 3.1 x 1O 9 M (1.0% of the missing baryons from the Galaxy), or a 
truncated NFW halo of mass 4.1 x 1O 9 M (1.4% of the missing baryons). Using our assumed 
metallicity of 0.2Z & , we recalculated the emissivity using the XSPEC utility with an APEC 
model (Smith et al. 2001) and all other parameters matched to the observations in Kuntz 
and Snowden (2000), which makes these halo masses increase by 190%. The corrected halo 
masses match very closely with the previous constraints: 5.9 x 10 9 M Q (2.0% of the missing 
baryons) for an NFW halo, 7.8 x 1O 9 M (2.6% of the missing baryons) for a truncated NFW 
halo. This halo has a cooling time larger than 10 Gyr for r > 15 kpc (no truncation) or 
r > 18 kpc (with truncation). 

We summarize all of these constraints in Table 1. Upper limits on the total hot halo 
mass, electron density at a radius of 50 kpc, and fraction of the missing baryons from the 
Galaxy in the halo, are computed for each constraint. The ranges on some parameters reflect 
ranges in the original constraint combined with fits to both a truncated and a non-truncated 
NFW profile. We also present the upper limits in terms of a flattened, power-law profile, 
discussed further in section 7. 
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Finally, we also include a lower limit on the Galactic hot halo in addition to the upper 
limits previously noted. Blitz and Robishaw (2000) measure HI mass for the dwarf spheroidal 
galaxies in the Local Group, and use these masses to place upper limits on the amount of 
ram-pressure stripping that these galaxies could have experiences. This leads to an upper 
limit on the density of the Local Group intragroup medium, corresponding to n e > 2.5 x 10~ 5 
cm -3 . The Galactic hot halo must be denser then the ambient medium around it, so this 
figure also serves as a strict lower limit on the density of the hot halo. 

Our picture of the Milky Way hot halo is therefore a halo obeying an NFW profile, 
and containing several x 1O 9 M of material. It extends out to at least 50 kpc, although 
its density declines to the IGM density by a radius of 80 kpc or so. Such a halo is fully 
consistent with all of the constraints we have identified, but it does not contain more than 
a few percent of the missing baryons from the Galaxy. 

4. OVII Absorption in Other Galaxies 

The previous techniques will not work for constraining the haloes of most other galaxies, 
since we cannot measure the DM of distant extragalactic pulsars, and even if we could, 
we cannot reliably separate the dispersion from intergalactic electrons from the dispersion 
caused by the galaxy's halo of missing baryons. However, the integrated column density of 
free electrons from the hot haloes still will have a measurable effect in other galaxies through 
OVII absorption features along quasar sightlines. 

At the temperatures of the haloes around large galaxies (a few hundred eV), intergalactic 
Oxygen exists in the highly ionized OVII or OVIII states, of which the former can be detected 
through its 21.60 A Ka absorption line if superimposed in front of a bright background X- 
ray source. Zero-redshift OVII absorption has been detected several times, and attributed 
variously to the Local Group IGM (Nicastro et al. 2002, Rasmussen, Kahn, and Paerels 
2003) or to a hot halo around the Milky Way (Wang et al. 2005, Fang et al. 2006, Bregman 
and Lloyd-Davies 2007). In addition to the local detections, there are a few published claims 
of nonlocal OVII absorption lines (Mathur, Weinberg, and Chen 2003, Nicastro et al. 2005b, 
Buote et al. 2009), but these have yet to yield any significant constraints (see, e.g., Kaastra 
et al. 2006). 

One successful approach comes from a statistical argument by Fang et al. (2006). They 
studied 20 active galactic nucleus (AGN) sightlines without any detections of nonlocal OVII 
absorption and argued based on the number density of L* galaxies and of Local Group- 
sized galaxy groups that most of the missing baryons must be restricted to haloes around 
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the former; if these haloes extended to megaparsec scales, it would be statistically unlikely 
(~ 3.7a) not to have seen them in OVII absorption. 

We extend this result by examining observations along specific sightlines, and broadening 
the search to galaxies smaller than L*. We select critical sightlines with the best combination 
of redshift interval and high signal-to-noise observations discussed by Fang et al. (2006), and 
count the intervening galaxies along the sightline, focusing on the twelve AGN that fall within 
the coverage of the Sloan Digitial Sky Survey (SDSS) Seventh Data Release (Abazajian et 
al. 2009). For each of these twelve sightlines, we searched the SDSS archive for foreground 
galaxies whose virial radii overlapped the sightline. If such a galaxy is present, but no OVII 
absorption is detected at its redshift, then we can place an upper limit on the OVII column 
from the galaxy. 

We searched within a 4° x 4° region around each AGN for SDSS galaxies with photomet- 
ric redshifts less than the AGN's redshift. We also found a known Lyman-a absorption line 
at each of these photometric redshifts to help verify the photometric redshifts; see Table 2. 
The photometric redshift was used to compute the distance to each galaxy and the angular 
distance between each galaxy and the AGN sightline. The virial radius for each galaxy was 
estimated, and galaxies were selected whose virial radii were larger than their distance to 
the AGN sightline. We inferred the virial radius by computing the SDSS g-band galaxy 
luminosity and assuming a scaling relation between luminosity and virial radius, normalized 
to the Milky Way (L g = 1 x 1O 1O L , R viT = 250 kpc). We consider two forms for the scaling 
relation. The first, (scaling relation A), is based on empirical scaling between the baryon 
fraction (taken to be a proxy for luminosity) and the total mass of a galaxy (McGaugh 2000, 
Dai et al. 2010): f b oc L/M oc M a5 , so L oc M 3 / 2 and R oc L 2/9 . This seems to represent 
the data over many orders of magnitude in galaxies with disk components (McGaugh et 
al. 2005). For elliptical galaxies, the well-established "tilt" in the fundamental plane (e.g. 
Tortora et al. 2009) has the baryon fraction scaling in the opposite direction with mass: 
M/L oc L 1 ^ 4 , so L oc M 4 / 5 and R oc L 5 / 12 (scaling relation B). We classify each galaxy as 
early-type or late-type using its SDSS u — r color, with early-type galaxies having u — r> 2.2 
and late-type galaxies having u — r < 2.2 (Strateva et al. 2001, Mateus et al. 2006). For 
the former, we apply scaling relation B, and for the latter we apply scaling relation A to 
estimate the virial radius. 

Using this technique, out of our sample of twelve sightlines we find three with nearby 
galaxies whose virial radii encompass the sightline: 3C273 (z em = 0.158), Q1821+643 (z em = 
0.297), and Ton 1388 (z em = 0.177). For each intervening galaxy, the mass and the baryon 
fraction are estimated, as well as the expected mass in missing baryons. We assume these 
baryons are distributed in a halo around the galaxy, extending to the virial radius, and that 
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they obey an NFW density profile (using several values for the concentration parameter 
between c=9 and c—17, as in Klypin, Zhao, and Somerville 2002). The expected electron 
column contributed by each galaxy is computed by integrating along the AGN sightline as 
it passes through the halo. This gives us N e>NFW) the electron column expected from a hot 
halo in hydrostatic equilibrium around each galaxy. 

In Table 2, we compare this column to the minimum detectable electron column inferred 
from Fang et al. (2006). For each sightline, they provide the expected width of a 3cr detec- 
tion of OVII absorption, from which we compute the OVII column that would produce this 
line using N(0 VII) = 3.48 x 10 14 EW, where N is in cm -2 and EW in mA (Bregman and 
Lloyd-Davies 2007). The electron column is estimated from the Oxygen column by adapting 
the conversion derived by Bregman and Lloyd-Davies (2007): 



where / is the O VII ion fraction (we use / = 0.5) and assuming gas of Z = 0.2Z Q . 
This yields N ejth , the minimum (threshold) electron column that can be detected at 3er using 
OVII absorption. These assumptions are conservative. The fraction of Oxygen in the O 
VII state is close to unity for temperatures between 3.5 x 10 5 K and 1.5 x 10 6 K, which are 
the expected temperatures for hydrostatic haloes around most of our selected galaxies. This 
is probably a low estimate since hot intracluster medium gas typically has a metallicity of 
about Z = 0.2Z Q , even towards the virial radius of distant galaxy clusters (Maughan et al. 
2008, Anderson et al. 2009). For both / and Z, choosing less conservative values reduces 
N e fh, which would imply even stronger expected detections than those listed in Table 2. 

Table 2 shows that two of these galaxies should have been detected in existing observa- 
tions. The nondetection of these absorption lines implies real electron column must be much 
smaller than the ^ e ,NF\v value derived above. In light of our conservative assumptions about 
unknown parameters like / and Z, the most reasonable way to reduce N e ^pw below N e>t h is 
to bring down the mass of the haloes. The radius could also be increased to several times the 
virial radius, which would reduce the density but increase the path length, but this would 
also increase the number of galaxies whose haloes intersect with the AGN sightlines. For 
example, if R is doubled to twice the virial radius, we would expect 4x as many absorbers, 
yielding about 6 detections of 3cr or more. Reducing the halo radius below the virial radius 
increases the density, which worsens the pulsar dispersion measure inconsistencies in the 
previous section and leads to a problem of excess X-ray emission (see next section). Also, 
our strongest constraints come from two galaxies of 5-10% of L*, each within about 40 kpc 




- 11 - 



of the sightline. As long as these galaxies have haloes at least 40 kpc in extent, we can put a 
3cr upper limit on the mass in the halo of 70% of the missing baryons, so reducing the radius 
over a reasonable range will not have much of an effect. Thus the only remaining option is 
to lower the mass, below 70% of the missing baryons at 3a. 

This is a stronger conclusion than Fang et al. (2006) reached in their analysis. They 
were able to rule out haloes with columns greater than with radii of 1 Mpc at 3.7a, but 
they did not rule out such haloes with R ~ 100 kpc. Their statistical analysis predicted only 
0-1 galaxies to fall within 100 kpc of a sightline in their entire sample, so their method did 
not really constrain the parameters of haloes of this size. Indeed, we find only one galaxy, 
and it lies at an impact parameter of 135 kpc, in agreement with their prediction. However, 
there are several 0.1L* galaxies closer than 100 kpc, and these are critical in deriving our 
constraints. At larger radii, our predictions are also consistent with not observing excess 
OVII absorption in a sightline that passes within 350 kpc of M31 (Bregman and Lloyd- 
Davies 2007). 



5. X-Ray Surface Brightness 

The failure to detect extended hot gaseous halos from X-ray emission (e.g. Strickland et 
al. 2004, Li et al. 2006, Li, Wang, and Hameed 2007) places particularly strong constraints 
on the fraction of missing baryons present around galaxies. Unlike the other constraints in 
our paper, this constraint has been known for at least a decade; see, for example, Benson 
et al. (2000) and Wu, Fabian, and Nulsen (2001). The former paper worked out X-ray 
luminosities of hot haloes around massive late-type galaxies, and found these should be 
some of the most luminous X-ray sources in the local Universe. The latter paper explored the 
unresolved X-ray background, and noted that the predicted contribution to the background 
from hot haloes around galaxy groups was an order of magnitude larger than the unresolved 
background. Since these haloes are so large and diffuse, however, measurements of both the 
luminosity and the unresolved background number density require extrapolation out to large 
radii and low surface brightness. Therefore, we instead consider this discrepancy in terms of 
X-ray surface brightness, since this quantity is directly measurable and can easily be used 
to establish observational constraints on these hot haloes. 

We calculate the X-ray emissivity for three different galaxy masses: 3 x 10 11 M Q , 8 x 
lO n M , and 2 x 1O 12 M . For each mass, the hot halo is assumed to obey an NFW profile 
with a concentration C=12 (the profile is not very sensitive to this value), and the hydro- 
static equilibrium temperature is computed. The emission from the baryons in the hot halo 
is modeled with an optically thin emission model (the APEC model) with photoelectric ab- 
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sorption from Galactic neutral Hydrogen of column density N# = 3 x 10 20 cm -2 . For each 
halo, we use the WebSpec utility to compute the flux F 1 over an energy range of 0.1-6 keV 
on the ACIS-S instrument. The metallicity was set to 0.2Z Q . 

We estimate a typical Chandra ACIS background of 1.3 x 1CT 12 erg s _1 arcmin~ 2 from, 
e.g., Bregman et al. (2009), which is converted into a flux (in erg s _1 cur 2 ) for each 
halo by dividing by the telescope effective area A determined using XSPEC. Assuming halo 
emission can be detected once it reaches a surface brightness of about 30% of the background, 
this background can be combined with the flux Fi from above and with the definition of 
normalization in the APEC model (at z — 0, norm = n 2 R x 10 _14 /47r) to calculate the 
minimum detectable emission measure for gas in each halo: 

, ~. An x 0.3 x (1.3 x 10~ 12 erg s _1 arcmin -2 ) - 

(n 2 R) min = r = — - — cm" 5 

10" 14 x Fl x (8.5 x 10- 8 arcmin sterad" 1 ) x A 

The results of this analysis are presented in Table 3. 

We integrate along columns through the halo to find the expected EM at each projected 
radius from this galaxy. In Figure 1, we compare the expected emission measure to the 
minimum detectable EM for galaxies of various concentrations at all three masses. While 
the emission within the disk may be obscured by or confused with the disk, the haloes can be 
detected out to many times the radius of the thick disk in all three cases. A Milky Way-sized 
halo, for example, should be detectable in emission out to about 75 kpc, and even a hot halo 
around a galaxy of total mass 3 x 1O 11 M should be detectable out to twice the radius of 
the thick disk. 

These haloes should be very obvious in X-rays around all reasonably-sized galaxies, and 
yet the observed X-ray haloes are confined to the immediate vicinity of the thick disk. This is 
a serious problem for the paradigm of hot haloes full of missing baryons. The missing baryon 
electron density must be smaller than the model predicts; the upper limit on the mass of the 
halo to bring the emissivity at the edge of the thick disk down to the detection threshhold is 
between 11% and 24% for the three different masses. In general larger galaxies have stricter 
upper limits, since they have more mass in their haloes. Therefore, this constraint can be 
improved further using deep X-ray observations of very massive isolated disk galaxies. 

Observations can also constrain the density profile of the hot halo; based on the predicted 
EM, in the detectable region outside the thick disk our NFW profiles roughly obey a (5 model 
with (3 ~ 0.87. This is slightly steeper than the profiles of hot gas in early-type galaxies, 
which have (5 ~ 0.4 — 0.6 (e.g. Brown and Bregman 2001). In Figure 5, we also consider 
a flatter density profile as discussed further in section 7. This profile has j3 ~ 0.36 in the 
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observable region, shallower than observed gas profiles. However, since observed values of 
(3 fall in between these two profiles, it seems likely that real haloes are neither as steep as 
NFW profiles nor as shallow as the flattened profiles of section 7. 

6. Energetics for Gas-Rich Late-Type Galaxies 

Having established some observational constraints on hot haloes around the Milky Way 
and other large galaxies, we comment briefly on the overall energetics of expelling the missing 
baryons from galaxies (either into virialized haloes or outside of the entire system) using 
large-scale galactic winds. These winds are a common component of the current generation 
of semi-analytic models of galaxy formation, required to reconcile the dominant "cold flows" 
picture (e.g. Keres et al. 2005), which assumes galaxies accrete the cosmological fraction of 
baryons, with the observed severe baryon depletion in nearly every galaxy today. Winds are 
also observed directly at z ~ 2 — 3 in galaxies undergoing periods of intense star formation 
and AGN activity (e.g. Pettini et al. 2001, Geach et al. 2005, Law et al. 2007), but appear 
insufficient to expel the missing baryons (Martin 2005). 

Such winds come in two varieties: "energy-driven" and "momentum-driven". The for- 
mer carry kinetic energy proportional to the energy injected into the interstellar medium by 
supernovae and/or an AGN, while the latter (Murray, Quataert, and Thompson 2005) carry 
momentum proportional to the momentum injected into the ISM by supernova shocks and 
radiation pressure from stars and AGN. In both cases, however, the constant of proportion- 
ality (the "efficiency" of the coupling) has to be large in order to drive enough of a wind to 
expel the missing baryons. For example, recent simulations by Dutton and van den Bosch 
(2009) require 25% efficiency for energy-driven winds or 100% efficiency for momentum- 
driven winds to expel the missing baryons, and they neglect cooling of the ejected baryons, 
which can result in the baryons re-accreting onto the disk in a galactic "fountain" and can 
therefore increase the required ejection energy by several times. Oppenheimer and Dave 
(2008) use momentum-driven winds with an equivalent efficiency of 200% with respect to 
the energy available to energy-driven winds. Such high efficiencies are an attractive feature 
of momentum-driven winds, but we note that the theory still needs to be fully developed and 
observationally supported. In particular, it is not yet clear that momentum-driven winds 
are stable enough to drive most of the baryons from a galaxy instead of fragmenting into a 
less efficient multiphase wind, although further research may illuminate this question. 

Regardless, while requiring such a large efficiency may already be problematic, driving 
these winds also leads to some challenging conclusions. There are two major sources of 
energy /momentum available for this purpose: supernovae and accretion onto a supermassive 
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black hole (SMBH). In addition, one or both of these can be triggered by major mergers, but 
mergers primarily increase the amount of energy in these two reservoirs rather than adding 
a new component. (Mergers also do release gravitational energy, but we are not aware of 
any way to couple efficiently this energy to the baryons in galaxies.) The maximum energy 
or momentum available in a galaxy from supernovae should be proportional to the stellar 
mass of the galaxy, and the maximum energy/momentum available from the AGN should be 
proportional to the mass of the black hole. Since the black hole mass is proportional to the 
galactic bulge mass (Marconi and Hunt 2003) , which is dominated by stars, this would seem 
to suggest that the total energy available to eject baryons is related to the stellar mass of 
the galaxy, which would lead naturally to a baryonic Tully-Fisher relation (McGaugh 2005). 

However, this picture also leads to the prediction that variations in the energy sources 
should affect the baryon fractions of galaxies, and this is not observed. The observations 
show that, at least for disk galaxies, the baryon fraction seems to depend almost solely on the 
circular velocity at large radius (a proxy for total mass; McGaugh 2005, Stark, McGaugh, 
and Swaters 2009). But if the baryon fractions are caused by supernova- and AGN-driven 
winds, then galaxies with fewer supernovae or smaller SMBHs should have higher baryon 
fractions than galaxies with more supernovae or larger SMBHs. So, for example, galaxies 
with no central black hole such as M33 (Gebhardt et al. 2001) should have higher baryon 
fractions than galaxies with the the same total mass but a large SMBH, when in reality M33 
lies on the BTF. 

If the supermassive black hole (SMBH) in a galaxy is responsible for a significant amount 
of the baryon loss, there should be a relationship between the baryon fraction and the galaxy 
type. Equivalently, in a fit to the baryon Tully-Fisher relationship, galaxies with relatively 
larger SMBHs (earlier-type galaxies) would lie below the normal scaling relation if AGN 
activity from a SMBH expels gas. To search for this effect, we examine the sample of galaxies 
from McGaugh (2005) and use the galaxy type as a proxy for the relative importance of a 
SMBH. When the galaxy type was not determined to sufficient precision, the galaxy was 
not used for this comparison. For the 60 galaxies with good galaxy types, there is no trend 
of the residual with galaxy type (Figure 2). We estimate that we could have detected a 
trend of about 20% or larger in the residual between Sb and Sd, suggesting AGN activity is 
responsible for removing a mass equivalent to less than 20% of the remaining baryons. Since 
present-day baryon fractions are typically 5-30% of the cosmological value, this corresponds 
to AGN removing no more than 1-5% of the original baryonic content from galaxies over 
their lifetime. 

We apply a similar analysis to examine the effect of supernova energy on baryon fraction. 
We computed the ratio of stellar mass to total baryon mass in the sample of gas-rich disk 
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galaxies of Stark, McGaugh, and Swaters (2009). To estimate the total baryon mass, we 
assumed an NFW halo potential and applied the relation M 2 oo = 2.3 x 1O 5 i>|/i _1 M for 
Vf in km s _1 from Navarro (1998). The stellar mass is computed in Stark et al. for each 
galaxy from stellar population modeling. In Figure 3, we compare this ratio to the circular 
velocity at large radius Vf for each galaxy. The figure shows these galaxies are all significantly 
depleted in stellar fraction compared to average disk galaxies, so there should be significantly 
less supernova energy available to remove baryons in these galaxies. We examined whether 
the residuals about the BTF are correlated with the degree to which the galaxies are gas- 
dominated. No such correlation was found, using both parametric and non-parametric tests, 
which suggests that supernovae from stars in these galaxies do not appear to be responsible 
for removing most of the baryons. 

An alternative explanation of the lack of BTF correlation with galaxy type or gas 
fraction can be inferred from the result of Navarro and Steinmetz (2000). They noted that 
the remarkably low scatter in the I-band Tully-Fisher relation can be explained through 
galaxies with different disk mass fractions scattering along the Tully-Fisher relation instead 
of perpendicular to the curve. The idea is that the disk mass actually contributes to the 
measured rotation speed in addition to the halo mass, so if the disk mass increases, the 
rotation speed will also increase, regardless of the halo mass. This explanation does not apply 
directly to our argument, however. Neither the gas fraction nor the galaxy type correlates 
with the disk mass, so there is no reason to expect the scatter between these quantities 
and the rotation speed to follow the BTF. However, if supernovae or AGN accretion has 
expelled a dynamically significant mass of baryons from one of these galaxies, we would 
expect the rotation speed of the galaxy to decrease near the disk, in an analogous effect to 
the result of Navarro and Steinmetz. To minimize this concern, the rotation speeds used 
in our analysis are inferred from 21 cm emission at large radii, outside of the disk, where 
the rotation curve has flattened (for details, see McGaugh 2005 and Stark et al. 2009). At 
these radii (tens of kpc), the dark matter dominates the gravitational potential, and so the 
efficiency of expelling baryons should not affect the measured rotation velocity. Thus we 
do not believe correlations between baryon fraction and rotation velocity affect our claim 
that AGN accretion and supernovae are not primarily responsible for expelling the missing 
baryons. 

A simple way to understand the argument of this section is to consider two galaxies: 
NGC 4183 and NGC 4138. Figure 4 shows images of these two galaxies taken from the 
Sloan Digital Sky Survey, and in Table 4 we present measured or inferred parameters for the 
baryons in these galaxies, based on the results in NED and in Stark, McGaugh, and Swaters 
(2009) as described above. 
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These two galaxies are similar in mass, but very different in composition. NGC 4138 is a 
bulge-dominated lenticular galaxy with 91% of its baryons in stars. NGC 4183 is a late-type 
disk galaxy with only 31% of its baryons in stars (the other 69% are in gas). Our sample 
contains more extreme cases of both gas-rich and gas-poor galaxies than these, but these 
two present a useful contrast in both gas fraction and galaxy type, with little difference in 
total mass. Since NGC 4138 has such a high stellar fraction and such a prominent bulge, 
it must have been able to use its AGN and its supernovae to produce a significant galactic 
wind and drive out more of its baryons than NGC 4183, and yet the latter has a smaller 
baryon fraction than the former. 

More quantitatively, if we assume a negligible SMBH in NGC 4183 and assume its 
2.1 x 1O 9 M of stars has produced 1.8 x 10 7 supernovae^, we can estimate the average mass 
each supernova needs to eject to account for the missing baryons. From eq. (1), we calculate 
NGC 4183 is missing 7.9 x 1O 1O M of baryons, so each supernova needs to heat and eject 
44OOM from the galaxy. It already requires 180% of the supernova energy (i.e., 1.8 x 10 51 
erg) to heat 4400M Q to the virial temperature of NGC 4183, and a comparable amount of 
energy to drive this material up into a halo from the gravitational potential of the disk, 
for a total supernova energetic efficiency of 350% (or more, if cooling is included). These 
extreme efficiencies may pose a problem even for momentum-driven winds, and this leads to 
our suggestion that the classical picture of AGN- and supernova-driven winds is incomplete 
for explaining the missing baryon problem. 



There exists one significant modification to the standard picture of hot haloes that must 
be considered. Significant pre-heating of the baryons can lead to excess entropy in the halo 
gas, which causes it to deviate from the dark matter NF W profile and obey a flatter density 
profile instead (e.g. Mo and Mao 2002, Kaufmann et al. 2009). We approximate the relaxed 
high-entropy profile in Figure 6 of Kaufmann et al. (2009) with a power-law: 



with p = 0.9. We recomputed the Galactic constraints for this flattened density profile (see 



7. The Effect of Flatter Density Profiles 




(3) 



1 This is estimated by integrating over a Chabrier (2003) initial mass function, letting stars with M > 8M Q 
explode in a core-collapse supernova, and applying a 20% correction (Mannucci et al. 2005) to include type 
la supernovae 
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Table 1). The DM constraint is significantly relaxed: the halo can hold 58% of the missing 
baryons before reaching the LMC DM limit of 23 cm -3 pc. This limit is dramatically weaker 
compared to the NFW upper limit of about 5%. On the other hand, the pressure confinement 
constraints based on high velocity clouds and Magellanic Stream clouds place an upper limit 
of 6 — 13% of the missing baryons (2 — 4 x 10 10 M Q ), which still fails to resolve the second 
missing baryon problem. A flattened profile containing all the missing baryons also violates 
the EM measurement of Kuntz and Snowden (2000). They measured an EM of 2.1 x 10~ 3 
cm -6 pc (which becomes approximately 7.7 x 10~ 3 cm" 6 pc for our assumed metallicity), 
while the flattened halo would have an EM of 1.9 x 10~ 2 cm" 6 pc. Thus the picture we 
described in Section 3.2 remains - the Galaxy has a limited hot halo of some form that 
declines to the IGM density, but it does not contain most of the missing baryons from the 
Galaxy. 

The extragalactic emission measure constraint (X-ray surface brightness) is also signif- 
icantly affected by a flatter profile. Kaufmann et al. claim their profile lowers the X-ray 
luminosity of the inner region of the halo by a factor of ~ 70. We computed the emission 
measure for the flattened power-law profile (Figure 5) and indeed found the EM lowered by 
a factor of 70-300 at the edge of the thick disk. However, the emissivity of the gas per unit 
mass increases, by about 75%, because the flattened profile is hotter than the NFW profile. 
Finally, since the slope of the profile is flatter, the emission can remain significant out to 
larger radii. The result of these combined effects is that the emission from within the disk 
is lowered significantly, but the emission remains flatter, so at larger radii it is still visible. 
As Figure 5 shows, smaller galaxy haloes are not detected, but a Milky Way-sized galaxy 
would still have a halo detectable out to 70 kpc or so, and a larger galaxy would have a hot 
halo visible to even larger radii. Searching for X-ray emission at significant fractions of the 
virial radius around large, isolated disk galaxies is therefore a good way to constrain these 
flattened profiles. 

This claim differs from the conclusions in Kauffmann et al. (2009). They argue that the 
surface brightness of their high-entropy halo is consistent with observations, while we argue 
the emission from such a halo should be detectable tens of kpc away from the galaxy disk. 
This distinction is primarily due to different assumptions about the baryon fraction. They 
initialize their models 10 Gyr ago with a baryon fraction f b = 0.1, while we assume galaxies 
form with the cosmological baryon fraction /& = 0.17. This gives our halo a density 1.7 
times larger than theirs, which yields a factor of three in the emission measure. Correcting 
for this effect reconciles our respective predictions with each other, although we disagree on 
the interpretation. To understand the second missing baryon problem, however, one must 
begin with the cosmological baryon fraction, and then explain the fate of all the baryons 
originally associated with the galaxy over cosmic time; thus we argue that flattened density 
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profiles do not hide all of the missing baryons from galaxies. 

The flattened profiles are still detectable through their OVII absorption as well. At 
large fractions of the virial radius, a flattened profile has a higher density compared to an 
NFW profile, so sightlines that only pass through a small section of the hot halo produce 
larger absorption lines if the profile is flatter. Sightlines that pass closer to the central galaxy 
produce roughly equivalent sightlines to the NFW halo. We computed the absorption lines 
expected for our seven absorbers, and displayed the results in Table 2 above. In the former 
regime, the small (< la) detections predicted for NFW haloes increased moderately if the 
profile was flatter and the 1 — 2a detections became roughly 3a detections. In the latter 
regime, the two large > 4a detections both decreased slightly. Since a 3a absorption line 
could be observable, the flattened density profiles do not seem to solve the OVII absorption 
issue, and actually make the problem worse. 

In general, flattened profiles can contain a larger fraction of the missing baryons from 
galaxies than NFW profiles. For the Milky Way, a flattened profile can meet the LMC pulsar 
DM constraint with up to 58% of the missing baryons in the hot halo instead of 5%, although 
the pressure constraints bring this figure down to 6-13%. The X-ray emissivity argument 
constrains the mass of these haloes around other large galaxies to about 59% of the missing 
baryons from these galaxies. This is a factor of three more mass than the NFW profile can 
hold without violating surface brightness limits. The OVII argument also limits the mass 
of these haloes in the same way it limits the mass of NFW haloes, with an independent 
upper limit of about 70% on the fraction of the missing baryons that can reside in the halo. 
Flattened haloes are therefore an improvement, but they do not seem to be a perfect solution. 

If the halo is driven from the galactic disk by a superwind, then the energetics problem 
we highlighted still applies. The flattened profile puts some additional constraints on the 
wind, however. In Kaufmann et al. (2009), the high entropy that keeps the profile flat is 
assumed to be injected at redshift 2 or 3. This pushes back the energetics problem 10 Gyr 
or so, but still does not quite resolve it. On the other hand, if the halo has existed in its 
present form for 10 Gyr, perhaps the gas simply never fell into the disk in the first place. 
We will discuss this idea further in the conclusion. 



8. Discussion and Conclusion 

The above arguments place separate and independent constraints on the baryons around 
galaxies. These arguments lead us to two general conclusions: 

First, the majority of the missing baryons from spiral galaxies do not lie in hot haloes 
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around these galaxies. For the Milky Way, the combination of observed constraints suggests 
that the common assumption of a hot halo of missing baryons obeying an NFW density profile 
with a concentration of 12 can only hold 2-3% of the baryons missing from the Galaxy. The 
emission measure of theoretical haloes places an upper limit of 11-24% on the fraction of the 
missing baryons that can reside in similar hot haloes around other large galaxies. 

The analysis of O VII sightlines constrains hot haloes around a wider mass range of 
galaxies. The free electrons in these haloes should be detected through O VII absorption 
in the spectra of several strong background X-ray sources. This absorption has not been 
observed, so the column must be several times smaller than the prediction. This would 
occur if the halo has a radius of < 40 kpc or if its mass is reduced by factor of 3-4. The 
former case leads to further difficulties. The density would be 0.01-0.1 cm -3 , so the halo 
would cool quickly and should be far too bright in X-rays. The latter case is possible, but 
such a rarefied halo would not account for the missing baryons from galaxies. 

The proposal of hot haloes obeying flattened density profiles instead of NFW profiles 
was also considered. Around the Milky Way, a typical flattened profile could hold 6 — 13% 
of the missing baryons without violating any of our constraints. The emission measure of 
hot haloes with a flattened profile around other large disk galaxies sets an upper limit on 
the mass of about 59% of the missing baryons from these galaxies. The OVII constraint is 
approximately unchanged, which also limits the mass and extent of flattened density profiles. 

Collectively, these arguments rule out the existence of a large reservoir of hot gas around 
large galaxies that could account for the missing baryons. The most likely scenario is the 
presence of a small gaseous halo around L* galaxies containing 10 9 — 1O 1O M and extending 
to a radius of 50-100 kpc. This picture also satisfies pulsar DM observations, measurements 
of the pressure around high- velocity clouds, and X-ray surface brightness constraints. If such 
a limited halo exists, it is not the primary reservoir of the missing baryons from galaxies. 

Our second general conclusion is that galaxies do not expel their baryons primarily 
by galactic winds driven by supernovae or AGN activity. We infer this from the lack of 
correlation between missing baryon fraction in galaxies and the stellar fraction or bulge mass 
in the galaxy. Since the latter observables are taken to be proxies for the energy available in 
a galaxy from supernovae or from a supermassive black hole, the lack of correlation implies 
these energy sources do not explain the observed missing baryon fractions. 

Additionally, there are galaxies where most of the baryons are still gaseous yet they 
have the same missing baryon fraction as their counterparts where the baryons are nearly all 
in stars. To explain the baryon depletion of some of these star-poor, gas-rich galaxies, the 
energy required is several times that expected from galactic winds powered by supernovae. 
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It is beyond the scope of this paper to model the remaining mechanisms for resolving 
the second missing baryon problem, but we can comment on the possible solutions. Based on 
our conclusions, a straightforward way to account for the missing baryons from galaxies is if 
the baryons never fell into galaxies in the first place. This would require pre-heating before 
or just as the galaxy is forming. The most obvious source of non-gravitational energy at 
these times is early supernovae, whose contribution to heating would be much more effective 
in the epoch of protogalaxies (when potential wells are shallow) than it is today. These 
early supernovae could even come from the hypothetical "Population III" stars, which have 
the advantage of a top-heavy IMF, so this process would not convert much of the missing 
baryonic mass into low-mass stars. 

If the missing baryons never accreted onto protogalaxies because of early supernova 
heating, there are several observable consequences in addition to resolution of the second 
missing baryon problem. These early supernovae pre-enrich gas in addition to pre-heating 
it, and this helps to explain phenomena like the G-dwarf problem and the "floor" in iron 
abundance noted in section 4. For more discussion of this possibility, see Benson and Madau 
(2003). Due to the top-heavy IMF, they would have short lifetimes, so few of these early 
unbound stars would survive today, and the strict limits on intracluster light (e.g Krick and 
Bernstein 2007) would not be violated. A top-heavy IMF would also increase the efficiency 
per unit mass of metal production, which has further implications for the metallicity in galaxy 
clusters, as discussed in Bregman, Anderson, and Dai (2010). There are also observable 
constraints on this process. For example, if supernova heating is to prevent the baryons 
from accreting onto galaxies, the supernovae must occur before the galaxies form, i.e. z > 6. 
There is also an important period at z ~ 8 when Compton cooling can efficiently cool the 
baryons (Cen 2003). Depending on the details of the supernovae and the cosmology, the 
bulk of the early star formation and baryon expulsion might occur either before or after 
this period. More detailed observation and theoretical work can further constrain these 
possibilities. 

The remarkably tight correlation for galaxies between dark matter halo mass and baryon 
fraction seems to offer a powerful clue for understanding galaxy formation. The correlation 
persists over variations in star/gas mass ratio The mechanism and details of this connection 
remain unclear, however, and still need to be understood. Further observations could also 
help to constrain more tightly the mass and extent of hot haloes around galaxies. Direct 
imaging around large galaxies in soft X-rays could strengthen the surface brightness con- 
straints. Measurements of extragalactic pulsars or additional OVII lines could also provide 
useful constraints on the electron column density around galaxies. Detailed examination of 
simulations could also be helpful, predicting, for example, the fraction of dark matter in 
the Universe that falls into galactic haloes for comparison to the behavior of the baryons. 



Finally, more modeling of galaxy formation could greatly assist in producing testable pre- 
dictions that would differentiate between various solutions to this second missing baryon 
problem. 
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Fig. 1. — Integrated X-Ray emission measures for a halo of hot gas obeying an NFW profile in 
hydrostatic equilibrium. The three colors denote the total mass of the host galaxy and its dark 
matter halo; each line represents the EM for an NFW profile of concentration parameter C=12. 
The three thick horizontal bars represent the approximate radial extent of the thick disk of a galaxy 
of the given mass; within these radii, the bright emission should be obscured and unobservable. 
The three thin horizontal lines are estimates of the minimum emission measures detectable with 
Chandra ACIS-S for each halo. These haloes are all observable out to many times the radius of 
the thick disk, and should be detected unless the halo mass is reduced to 11-24% of the mass of 
missing baryons. The slope of the surface brightness between the thick disk and the detection limit 
follows a /3-model with f3 ~ 0.87 in each case. 
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Fig. 2. — Comparison of the scatter of galaxies around the baryonic Tully-Fisher relation (BTF), 
based on the sample from McGaugh (2005). We use the galaxy type as a proxy for the size of the 
supermassive black hole (SMBH). For the 60 galaxies with good galaxy types, there is no trend of 
the residual around the BTF with galaxy type. As discussed in section 6, we conclude that, at least 
for galaxies of type Sb and later, AGN activity associated with SMBHs does not play a dominant 
role in the removal of baryons from those galaxies. 
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Fig. 3. — Comparison of galaxy residuals about the baryonic Tully-Fisher relation as a function of 
their stellar mass fraction for the gas-dominated sample of Stark et al. (2009). For several of these 
galaxies, the stars comprise less than 25% of the baryons, yet they have the same baryon fraction as 
the more star-rich galaxies. As discussed in section 6, we conclude that supernovae from the stars 
in these galaxies are not responsible for ejecting most of the missing baryons from the galaxies. 
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Fig. 4.— SDSS images of NGC 4138 (left) and NGC 4183 (right). NGC 4138 is at a distance of 
15.4 ± 2.3 Mpc, NGC 4183 is at a distance of 17.3 ± 0.9 Mpc. Both images have a scale of 0.3961 
"/pixel. For more information on and discussion of these galaxies, please see section 6. 



Table 1: Constraints on the Milky Way Hot Halo 
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Parameters determined from fitting profiles to the constraints noted in section 3.2. frac is the ratio of the mass of the hot halo to the mass of the missing baryons 
from the Galaxy (3 X lO 1 ^©). 



Table 2: QVH Sightline Analysis 
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The impact parameter b is the distance between the galaxy and the relevant AGN sightline. Each galaxy was matched to Lyman a absorption detected at or very 
near its photometric redshift; the equivalent with is noted. The neutral Hydrogen column densities implied by the Ly-a lines arc ~ 10 13 - 10 14 cm - significantly 
lower than the electron column densities, as would be expected from highly ionized plasma. The Lyman a references arc: (a) Scmbach et al. 2001, (b) Sembach 
ct al. 2004, (c) Bahcall et al. 1993. The u — r color is used to distinguish between early-type (u — r > 2.2) and late-type (u — r < 2.2) galaxies when computing 
N e ,NFW- EWih is the equivalent width of a 3<r OVII detection for the Chandra observation of this AGN in Fang et al. (2006), and N e t ^ is the electron column 
corresponding to this EW. N ej jvFW is the electron column we expect from the hot halo around this galaxy if it contains the missing baryons (the uncertainty 
reflects la uncertainties propagated from the redshift, an assumed 5 kpc uncertainty in the impact parameter, the g-band photometry, scaling relation, and different 
choices for C in the NFW profile). (S/N)nf\v is the strength of the detection we would expect for this halo if it extended to the virial radius. N e p ; and SN p j are 
the column density and signal-to-noisc of a flatter, power-law profile that will be discussed in section 7. 
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Table 3: Halo Emissivity Analysis 
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Numbers used for derivation of emission measure (n 2 i?) m i n of halo gas in hydrostatic equilibrium that would be detected by 
Chandra ACIS-S at 30% of background. Mass is the total mass of the galaxy and dark matter halo, and kT is the temperature 
of the halo gas. -Fi is the estimated flux from an APEC model (see text). And A is the effective area of the detector for each 
observation. 



Table 4: Comparison of Two Nearby Disk Galaxies 



Name 


i 
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Type 
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NGC 4138 
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9.30 


10.30 


SO 
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9.85 


0.016 


9.69 


9.33 
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Data taken from Stark, McGaugh, and Swatcrs (2009) and from the NASA Extragalactic Database, as described above. 
M200 is computed from Vf as described above as well, assuming h = 0.72. 
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Fig. 5. — As in Figure 1, the integrated X-Ray emission measures for a halo of hot gas, but in 
this plot the gas obeys a power-law profile: n(r) oc r~ ' 9 . These hot haloes have a lower emission 
measure, but a higher emissivity per gram due to their increased temperature. The largest should 
be detected unless the halo mass is reduced to 59% of the mass of missing baryons, the smaller 
haloes have weaker constraints. The slope of the surface brightness between the thick disk and the 
detection limit follows a /3-model with f3 ~ 0.36 in each case. 



